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Abstract 

Future colliders will probe the electroweak theory at energies much larger 
than the gauge boson masses. Large double (DL) and single (SL) logarithmic 
virtual electroweak Sudakov corrections lead to significant effects for observable 
cross sections. Recently, leading and subleading universal corrections for ex- 
ternal fermions and transverse gauge boson lines were resummed by employing 
the infrared evolution equation method. The results were confirmed at the DL 
level by explicit two loop calculations with the physical Standard Model (SM) 
fields. Also for longitudinal degrees of freedom the approach was utilized for 
DL-corrections via the Goldstone boson equivalence theorem. In all cases, the 
electroweak Sudakov logarithms exponentiate. In this paper we extend the same 
approach to both Yukawa enhanced as well as subleading Sudakov corrections 
to longitudinal gauge boson and Higgs production. We use virtual contributions 
to splitting functions of the appropriate Goldstone bosons in the high energy 
regime and find that all universal subleading terms exponentiate. The approach 
is verified by employing a non-Abelian version of Gribov's factorization theorem 
and by explicit comparison with existing one loop calculations. As a side result, 
we obtain also all top- Yukawa enhanced subleading logarithms for chiral fermion 
production at high energies to all orders. In all cases, the size of the subleading 
contributions at the two loop level is non-negligible in the context of precision 
measurements at future linear colliders. 
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1 Introduction 



The high energy behavior of the Standard Model (SM) will become increasingly im- 
portant at future colliders investigating the origin of electroweak symmetry breaking. 
At the expected level of precision required to disentangle new physics effects from the 
SM in the O (< 1%) regime, higher order electroweak radiative corrections cannot be 
ignored at energies in the TeV-range. The largest contribution is contained in elec- 
troweak double logarithms (DL) of the Sudakov type and a comprehensive treatment 
of those corrections is given in Ref. j l| to all orders. The effects of the mass-gap 
between the photon and Z-boson has been considered in recent publications || |3| since 
spontaneously broken gauge theories lead to the exchange of massive gauge bosons. 
In general one expects the SM to be in the unbroken phase at high energies. There 
are, however, some important differences of the electroweak theory with respect to an 
unbroken gauge theory. Since the physical cutoff of the massive gauge bosons is the 
weak scale M = Mw ~ Mz ~ Mh, pure virtual corrections lead to physical cross 
sections depending on the infrared "cutoff". Only the photon needs to be treated in 
a semi-inclusive way. Additional complications arise due to the mixing involved to 
make the mass eigenstates and the fact that at high energies, the longitudinal de- 
grees of freedom are not suppressed. Furthermore, since the asymptotic states are not 
group singlets, it is expected that fully inclusive cross sections contain Bloch-Nordsieck 
violating electroweak corrections [|]]. 

It has by now been established that the exponentiation of the electroweak Sudakov 
DL calculated in Ref. [p]] via the infrared evolution equation method with the fields 
of the unbroken phase is indeed reproduced by explicit two loop calculations with the 
physical SM fields [fj, [| ||. One also understands now the origin of previous disagree- 
ments. The results of Ref. Q , based on fully inclusive cross sections in the photon, is 
simply not gauge invariant as already pointed out in Ref. [|TJ. The factorization used 
in Ref. || is based on QCD and only takes into account contributions from ladder 
diagrams. In the electroweak theory, the three boson vertices, however, do not simply 
cancel the corresponding group factors of the crossed ladder diagrams (as is the case in 
QCD) and thus, infrared singular terms survive for left handed fermions (right handed 
ones are effectively Abelian) in the calculation of Ref. ||. The infrared evolution equa- 
tion method does not encounter any such problems since all contributing diagrams are 
automatically taken into account by determining the kernel of the equation in the 
effective regime above and below the weak scale M. It is then possible to calculate 
corrections in the effective high energy theory in each case yielding the same result as 
calculations in the physical basis. Thus, the mass gap between the Z-boson and the 
photon can be included in a natural way with proper matching conditions at the scale 
M. For longitudinally polarized gauge bosons it was shown in Ref. || that the DL 
kernel can be obtained from the Goldstone boson equivalence theorem. 

The picture that emerges has a clear physical interpretation. At high energies, 
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where particle masses can be neglected, the effective theory is given by an unbroken 
SU{2) x £7(1) theory for fermions and transversely polarized gauge bosons, and by the 
equivalence theorem for longitudinally polarized gauge bosons. The contribution from 
soft photons and collinear terms below the weak scale is determined by QED (including 
mass terms in the corresponding logarithms). This approach was utilized in Ref. |J to 
obtain the subleading (SL) universal terms to all orders for external fermion lines (up 
to Yukawa enhanced terms) and transversely polarized gauge bosons. Universal, i.e. 
process independent, are those terms which through Ward identities can be related 
to external lines and at high energies are given by the contributions to the virtual 
splitting functions (see Ref. ||). In addition, there are non- universal angular terms 
of the type log j log which can be important and should be included at least at 
one loop. In general, these terms cannot be resummed, they are process dependent 
and don't factorize with respect to the Born amplitude. At one loop, however, there 
is a general method for calculating such terms |]IU| and in practice at most a two loop 
approach to subleading logarithmic accuracy would be needed. 

An important aspect of resumming universal terms is given by the fact that there 
is a partial cancellation between the DL and SL corrections at energies around a TeV, 
thus enabling one to see how reliable a DL analysis for a given process really is. In 
addition, these are predictions of universal terms which can always by used to check 
higher order calculations which in the electroweak theory are extremely involved due to 
the number of mass terms and diagrams contributing. It is also conceptually important 
for a theoretical understanding of the infrared behavior of the SM. By comparing 
the subleading universal terms with existing one-loop calculations, we gain further 
evidence for the overall method employed, in particular when it comes to understanding 
differences between unbroken and broken gauge theories. 

From a phenomenological point of view, the corrections to longitudinal gauge bo- 
son production are important in case of a strongly interacting W sector without a 
fundamental Higgs boson. Our perturbative approach would of course break down in 
that case, however, it is important to know the precise form of the deviation of the 
new dynamics from the SM in the TeV range in order to understand the new physics 
behind the electroweak breaking sector. Corrections to Higgs bosons are important for 
a precise measurement of the Yukawa couplings at high energies in order to establish 
the Higgs mechanism. At the level of 6 — 8%, these corrections can certainly not be 
neglected at 1 TeV for determinations of the top- Yukawa coupling at e + e~ colliders 
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In this paper, we complete the all orders resummation of all SL universal Sudakov 
corrections to the SM. While in Ref. have restricted ourselves to calculating terms 
analogous to QCD, we now consider terms typical for broken gauge theories. These 
are in particular longitudinal degrees of freedom, processes with external Higgs bosons 
and Yukawa enhanced logarithmic corrections. 

In section El we discuss how SL contributions are obtained in the scalar sector and 



2 



the application of the equivalence theorem. In section [3] we give results for the vir- 
tual contributions to splitting functions involving Goldstone and Higgs bosons fulfilling 
evolution equations analogous to the Altarelli-Parisi equations. The correctness of this 
approach to the SL level is verified in section f| by employing a non-Abelian version 
of Gribov's bremsstrahlung theorem to processes involving Yukawa enhanced contri- 
butions. A similar approach can be utilized to verify analogous corrections for chiral 
fermions (6l,^l,^r) in section j|. Semi-inclusive cross sections for physical observables 
are given in section ^| and we compare our results with existing one loop calculations 
in section [7|. We discuss the size of the corrections obtained in section || and present 
our conclusions in section [|. 



2 The effective high energy theory 

In the following we discuss the corrections in a non-Abelian theory with scalar external 
"quarks", i.e. external scalar bosons charged under the unbroken gauge group. The 
physical picture is that at high energies, we can use this effective theory via the Gold- 
stone boson equivalence theorem to describe the longitudinal degrees of freedom. The 



latter will be discussed in detail in section 2.2 . An additional complication is given 
by the presence of Yukawa enhanced logarithmic corrections which are a novel feature 
of theories with spontaneously broken gauge symmetry. These terms are discussed in 
section |3|. We begin with a discussion of scalar (massless) QCD at high energies. 



2.1 Scalar QCD 

In this section we are interested in the collinear corrections to external scalars in a 
non-Abelian gauge theory at the subleading level when the external legs are taken on 
the mass shell. According to the discussion in Ref. [[J we have to calculate terms con- 
tributing to anomalous scaling violations. For purely virtual corrections the invariant 
matrix element fulfills the following differential equation for massless scalar quarks and 
all invariants 2pjpi ~ s large compared to an infrared cutoff [i and denoting t = log 



'| + ^ + fl ,(r. (t) -i^) + n.(r.( t ) + i7. 



ss 



xM(pi,...,p n ,g s ,v) = (1) 

to the order we are working here and where A4(pi, ...,p n , g s , fi) is taken on the mass- 
shell. The F(t) are infrared singular anomalous dimensions leading to DL corrections 
given in Ref. || while the gluon and scalar quark anomalous dimensions describe SL 
contributions. At higher orders, the subleading RG corrections can be incorporated 
by including a running coupling in each loop ]12]. An additional non-mass suppressed 



term occurs in the electroweak theory in the case of four scalar scattering amplitudes 
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Figure 1: A Feynman diagram determining the DL and SL contributions to scalar 
quarks in the on-shell scheme. In the massless theory there are scaling violations from 
loop corrections which can be described by anomalous dimensions. 



(such as — > due to the mass ratios in the coupling (A ~ -rjf) at the 

Born level. In that case one has to renormalize not only the gauge couplings but 
also the respective scalar couplings at the one loop level ||10|| . At higher orders, mass 
renormalization corrections are subsubleading. 

The universality of the corrections follows from the spin independence of the fac- 
torization theorems [13] and thus, from the universality of the splitting functions in 



scalar QCD[|. Such an equation is of course not physical for a theory with massless 
gauge bosons, however, in the electroweak theory purely virtual corrections can lead 
to physical cross sections. Only real soft photon emission needs to be included. 

We are interested in determining the scalar quark anomalous dimension 7 s3 which 
gives rise to logarithmic corrections due to scaling violations from the classical case. 
The factor of | in Eq. ([]]) originates from the fact that it is written for each external 
leg separately. 

Thus we must at one loop calculate the corrections depicted in Fig. [I] with the 
corresponding Feynman rules for scalar QCD: 

(1) b 3 r d n l (-Ihh + 2l(h - k 2 ) + Z 2 )(2Z + h- k 2 ) u 

v ~ F9s J (2n) n (I 2 - \ 2 + ie)((l + + ie)((1 - k 2 ) 2 + is) U 

This is in complete analogy to the situation in QCD |15, Requiring that the self 
energy corrections vanish on the mass shell and the on-shell vertex for zero momentum 



-"l would like to thank J. Collins for helpful discussions on this point. 
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transfer^ we find the following one loop result: 



M(1) = MBoni { 1+CF t* (4 log2 £ + 2 10 4 - 8 + 2 y) } (3) 

where A4 BoTn = —ig s T b (ki — k2) v . It is important to point out that both leading and 
subleading logarithmic corrections factorize with respect to the same group factor. The 
reproduction of this fact in the electroweak theory when compared to exact calculations 
with the physical fields is crucial in establishing the overall correctness of our approach 
since we must obtain the factorized form with respect to the effective high energy 
theory. The finite terms in Eq. are of course irrelevant to the discussion here and 
the infrared divergent soft and collinear terms were regulated using a fictitious gluon 



mass term. The difference to the QED result |T7| for the on-shell form factor 

aCo = "SB {i + £ H '° g2 £ + \ ' og ^ - 2 + 2 fl } < 4 > 

is (besides the coupling) mainly present in the different collinear divergent subleading 
term. This term differs due to the different spin of the particle emitting the gauge 
boson. Here .Mqed — — ie(A; 1 , T| r y^|^2, ' r ) as usual and replacing e 2 with Cpgl we 
obtain the QCD result from Eq. (^). Scaling violations for S-matrix elements can 
be described by calculating the anomalous dimension of the relevant gauge invariant 
operators. This is due to the fact that for massless theories there is a one to one 



correspondence between high and low energy scaling |L8], [OJ. Thus for the subleading 
scaling violations only the regions of large loop integration I in Eq. (|H) are relevant here 
(the double logarithms lead to infrared singular anomalous dimensions ptlfl ) and the 
corresponding anomalous dimension can be read off from the subleading logarithmic 
term: 

d ex 

Iss = TT, -{SaS + S a ) = ~C F - (5) 

O log /i 7T 

where 5 S s denotes the counterterm from the diagram depicted in Fig. [l], while S s 
corresponds to the wave function renormalization counterterm of the scalar quark. 
The sum in Eq. |5] is gauge independent. Since the factorization theorems of Refs. 
T3| do not depend on the spin of the quark, we can resum the leading and subleading 



virtual logarithmic corrections by using the Altarelli-Parisi equations. To this end we 
must formulate the above results in terms of the language of the splitting functions for 
a massless scalar quark. This will be done in the next section. First, however, we are 
going to discuss the scalar high energy sector in the Standard Model. In particular, 
there are additional corrections of the Yukawa type which need to be discussed, for 
which there is no analogue in unbroken gauge theories. 

2 We can perform this renormalization here since the non-Abelian components don't enter for the 
scalar quark anomalous dimension. The corresponding counterterm includes automatically the wave 
function renormalization contribution. 



5 



2.2 The equivalence theorem 

At high energies, the longitudinal polarization states can be described with the polar- 
ization vector 

e u L (k) = k u /M+0(M/E k ) (6) 

The connection between S-matrix elements and Goldstone bosons is provided by the 
equivalence theorem |]21|j . It states that at tree level for S-matrix elements for longi- 
tudinal bosons at the high energy limit M 2 /s — ► can be expressed through matrix 
elements involving their associated would be Goldstone bosons. We write schematically 
in case of a single gauge boson: 

M(W±^ phys ) = M^^^ + ot^) (7) 



M(Z L ,ip phys ) = iM{x,ip P hy S )+0 



The problem with this statement of the equivalence theorem is that it holds only at tree 
level p2| , f23|j . For calculations at higher orders, additional terms enter which change 
Eqs. (0) and (§). 

Because of the gauge invariance of the physical theory and the associated BRST 
invariance, a modified version of Eqs. (|7j) and (|3|) can be derived which reads 

k»M(W±(k),iP phys ) = C w M w M(<^(k),^) + o(^\ (9) 



k v M(Z u (k),iP phys ) = tC z M z M( X (k)^ phys ) + o(^j (10) 

where the multiplicative factors C w and C z depend only on wave function renormal- 
ization constants and mass counterterms. Thus, using the form of the longitudinal 
polarization vector of Eq. (|6]) we can write 

M(w£(k),*l) phya ) = C^M{4>Hk),^) + o[^ (11) 

M(Z L (k),^j phys ) = iC z M( X (k),^ phys ) + o(^\ (12) 

We see that in principle, there are logarithmic loop corrections to the tree level equiv- 
alence theorem. The important point in our approach, however, is that the correction 
coefficients are not functions of the energy variable s: 

C w = C w {-p, M, g, g>) , C z = C z (jt, M, g, gf) (13) 

The pictorial form of the Goldstone boson equivalence theorem is depicted in Fig. 
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Figure 2: The pictorial Goldstone boson equivalence theorem for W-pair production 
in e + e~ collisions. The correct DL-asymptotics for longitudinally polarized bosons are 
obtained by using the quantum numbers of the charged would be Goldstone scalars at 
high energies. 

|2] for longitudinal H^-boson production at a linear e + e~ collider. In the following we 
denote the logarithmic variable t = log where /x is a cutoff on the transverse part 
of the exchanged virtual momenta k of all involved particles, i.e. 

1? < k\ = mm(2(k Pl )(k Pj )/(p lPj )) (14) 

for all j / I. The non-renormalization group part of the evolution equation at high 
energies is given on the invariant matrix element level by ||: 

^- t M(L(k),% hys ) = K(t)M(L(k),^ phys ) (15) 

and thus, after inserting Eqs. (|1~T|), (p~2| ) we find that the same evolution equation also 
holds for 7V1 (0(A;), ^phys)- The notation here is L = {W^,Zl} and = {0 ± ,x} ) re_ 
spectively. Thus, the log -% dependence in our approach is unrelated to the corrections 
to the equivalence theorem, and in general, is unrelated to two point functions in a 
covariant gauge at high energies where masses can be neglected. This is a consequence 
of the physical on-shell renormalization scheme where the MS renormalization scale 
parameter JI ~ M. Physically, this result can be understood by interpreting the cor- 
rection terms C w and C z as corrections required by the gauge invariance of the theory 
in order to obtain the correct renormalization group asymptotics of the physical Stan- 
dard Model fields. Thus, their origin is not related to Sudakov corrections. In other 
words, the results from the previous section should be applicable to the subleading 
scalar sector in the electroweak theory regarding a non-Abelian scalar gauge theory as 
the effective description in this range. The only additional complication in the Stan- 
dard Model is the presence of subleading Yukawa enhanced logarithmic corrections 
which will be discussed below. It is also worth noticing, that at one loop, the authors 
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of Ref. [[It]] obtain the same result for the contributions from the terms of Eq. (p!3|). In 
their approach, where all mass-singular terms are identified and the renormalization 
scale ~p = y'i, these terms are canceled by additional corrections from mass and wave 
function counterterms. At higher orders it is then clear that corrections from two point 
functions are subsubleading in a covariant gauge. 



3 Subleading corrections from splitting functions 

In an axial gauge, collinear logarithms are related to corrections on a particular ex- 



ternal leg depending on the choice of the four vector n v [[14]]. In a covariant gauge, 
the sum over all possible insertions shown in Fig. |3] with all invariants large (~ s) 
is reduced to a sum over all n-external legs due to Ward identities. Overall, these 
corrections factorize with respect to the Born amplitude. We can therefore adopt the 
strategy to extract the gauge invariant contribution from the external line corrections 
on the invariant matrix element at the subleading level. In Ref. || we showed that 
in the high energy regime, subleading logarithmic corrections in massless theories are 
of collinear or RG origin. This is important since it allows us to use the Altarelli- 
Parisi approach to calculate the subleading contribution to the evolution kernel of Eq. 
(|j~5|). We are here only concerned with virtual corrections and use the universality 
of the splitting functions to calculate the subleading terms. For longitudinal degrees 
of freedom we have shown that to logarithmic accuracy the electroweak theory can 
be described by scalar Goldstone bosons via the equivalence theorem. Thus at high 
energies, the effective theory is analogous to scalar QCDQ. For this purpose we use 
the virtual gauge boson contributions to the splitting functions Pj±^,±(z), P^ x {z) and 
P BB {z) describing the probability to emit a soft and/or collinear virtual particle with 
energy fraction z of the original external line four momentum. The infinite momentum 
frame corresponds to the Sudakov parametrization with lightlike vectors. In general, 
the splitting functions Pba describe the probability of finding a particle B inside a 
particle A with fraction z of the longitudinal momentum of A with probability Pba to 



first order 15 



dV BA (z) = -^P BA dt (16) 

Z7T 



where the variable t = log -% for our purposes. It then follows 15 that 

fx 
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at, i \ oi s z{l-z)^r \V A —>b+c\ , 2 / nf7 N 
dV BA (z) = — ^ p d\ogk ± (17) 

Z7F Z spins K ± 



3 Although Yukawa terms are not present in QCD with scalar quarks, we will show in the next 
section that at subleading level the Yukawa terms can be treated as an additional term in the Altarelli- 
Parisi splitting function for Goldstone bosons. 
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where Va >b+c denotes the elementary vertices and 



z(l-z) 



E 

spins 



\V A - 



*B+C\ 



The upper bound on the integral over dk\ in Eq. (|T7| ) is s and it is thus directly related 
to dt. Regulating the virtual infrared divergences with the transverse momentum cutoff 
as described above, we find the virtual contributions to the splitting functions for 
external Goldstone and Higgs bosons: 



4,±4, 



Ti(Ti + 1) +tan 2 9 




-21og- 



2M 2 



6(1 -z) (19) 



The functions can be calculated directly from loop corrections to the elementary pro- 

25] and the logarithmic term corresponds to the 



16 



cesses in analogy to QCD [ 2l[ 
leading kernel of Ref. |§. We introduce virtual distribution functions which include 
only the effects of loop computations. These fulfill the Altarelli-Parisi equations^ 



dcf>(z,t) 
dt 



1 %{z/y,t)P^{y) 

y 



(20) 



The splitting functions are related by P^ = PE, + P^, where R denotes the contri- 
bution from real boson emission. P^ is free of logarithmic corrections and positive 
definite. The renormalizations with respect to the Born amplitude as well as the ones 
belonging to the next to leading terms at higher orders will be indicated below by 
writing a(p 2 ) — > a(s). 

Inserting the virtual probability of Eq. (|1^) into the Eq. ( pUP we find: 



0(1,*) 



) exp 



9_M 

8tt 2 



T,(T i + l)+tan 2 # w M log 2 4 





3 m? s 



(21) 



These functions describe the total contribution for the emission of virtual particles 
(i.e. z = 1), with all invariants large compared to the cutoff fi, to the densities <j>(z,t) 
{(f) = {0 ± , Xi H}). The normalization is not per line but on the level of the cross section. 



4 Note that off diagonal splitting functions do not contribute to the virtual probabilities to the 
order we are working here. In fact, for virtual corrections there is no need to introduce off-diagonal 
terms as the corrections factorize with respect to the Born amplitude. The normalization of Eq. 
( |l9[ ) corresponds to calculations in two to two processes on the cross section. The results, properly 
normalized, are process independent. 
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Figure 3: Feynman diagrams contributing to the infrared evolution equation (HD for 
a process with n external scalar quarks. In a general covariant gauge the virtual 
gluon with the smallest value of k± is attached to different external lines. The inner 
scattering amplitude is assumed to be on the mass shell. 



For the invariant matrix element involving external scalar particles we thus find at 
the subleading level: 



M(pi, ...,p n ,g s ,fi 2 ) = M(px, ...,p n ,g s {s))exp 




where 



9 2 {s) 



16tt 2 



Ti(Ti + 1) +tan 2 # 




4-41og4 



+ 



3 m\ 
2M2 



(22) 



log 



(23) 



Again we note that the running coupling notation in the Born-amplitude of Eq. fl23|) 
denotes the renormalization corrections of the Born amplitude and higher order cor- 
rections should be included by inserting a running coupling as in QCD The 



functions Wf correspond to the probability of emitting a virtual soft and/or collinear 
gauge boson from the particle subject to the infrared cutoff /i. Typical diagrams 
contributing to Eq. (|23|) in a covariant gauge are depicted in Fig. ^. The universality 
of the splitting functions is crucial in obtaining the above result. 

In addition to the Sudakov corrections in Eq. (p3|) we also have to include terms 
corresponding to the renormalization of the mass terms in the Yukawa coupling of 

the Born amplitude 



~ jfi, jfij at the one loop level jlQ| . At higher orders, mass 
renormalization terms are connected to two point functions and thus subsubleading. 
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4 Subleading corrections from Gribov's factoriza- 
tion theorem 



In this section we present further evidence for the results of the previous section by 
employing a non-Abelian version of Gribov's factorization theorem ||26|| . While in Ref. 



26] only real bremsstrahlung corrections are discussed, the form of the virtual soft 



and collinear divergences must factorize analogously due to the KLN-theorem [^7], |28 
The non-Abelian version is discussed in Ref. [|J. The essential point is that to DL 
accuracy the soft and collinear gauge boson with the smallest factorizes. From the 
definition in Eq. ( |i~4"|) it is clear that k\ fa \k\ 2 9 2 , where 9 denotes the angle between 
the emitted gauge boson of momentum k and the external line emitting this boson, i.e. 
both soft as well as collinear emission contributions are regularized simultaneously. In 
the Feynman gauge we then have 

% n r d^k p 'Pi 

M(p u ...,p n ;fi 2 ) = M Bo m(pi,-,Pn) - -77T\z Yl / o , i2 i • Ti \Ti \ 

2 (27r) 4 Js^-k 2 ^^ k 2 + it {kp j )(kpi) 

xT a (j)T a (l)M( Pl ,..., Pn ;k 2 ± ), (24) 

From Eq. ( |i~4"|) it is clear that ^^^^ = -j^- and that Eq. (|24] ) has the required factor- 
ized form. For the DL corrections it is convenient to employ the Sudakov parametriza- 



tion |29| given by: 

k = vpj + upi + k±_ (25) 



For the boson propagator we use the identity 



V -y + ir6(suv - k\) (26) 



k 2 + is suv — k ± + ie suv — 

writing it in form of the real and imaginary parts (the principle value is indicated by 
V). The latter does not contribute to the DL asymptotics and at higher orders gives 
subsubleading contributions. The cutoff will be introduced via the function Q(k 2 ± — fi 2 ). 
Rewriting the measure as d A k = d 2 k±d 2 k\\ with 

d 2 k ± = \k ± \d\k ± \d<p = ^dk 2 ± dip = irdk 2 ± (27) 
d 2 k\\ = \d(k° , k x ) /d(u,v)\dudv = \pj pi x — Pi Pj x \dudv ~ -dudv (28) 

where we turn the coordinate system such that the Pj,pi plane corresponds to 0, x and 
the y, z coordinates to the k± direction so that it is purely spacelike. The last equation 
follows from p 2 = 0, i.e. p 2 x ~ p 2 Q and 

(PjoPu ~ PhPj,) 2 ~ (PjoPk - PLPjJ 2 = iPiPif = ( s / 2 ) 2 (29) 
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Using in addition the conservation of the total group charge 



we arrive at: 



52T a {j)M(p h ..., Pj , ..., Pn , k 2 ± ) = (30) 
i=i 



2g 2 ™ f dk\ f 1 dv 



Zg f s 

M(p 1 ,...,p n ;n 2 ) = Mb ™(pi, -,Pn) ~ TT^Y, L , ■> I 

( 47r ) f^i -V fc_L J \*>±\/VS v 

x C l M(pi,...,p n ]k 2 ± ) , (31) 

where C\ is the eigenvalue of the Casimir operator T a (l)T a (l) (In QCD G\ = Ca for 
gauge bosons in the adjoint representation of the gauge group SU(N) and Ci = Cf 
for fermions in the fundamental representation). The differential form of the infrared 
evolution equation follows immediately from ([31]): 

dM{pi,...,p n ;fi 2 ) 2 2 

m 2 r = K(p)M(pi,...,p n ;fM), (32) 

where 

^ (/i) = "2g aiog(^) (33) 

with 

^W) = j^Q log 2 ± . (34) 

Wi is the probability to emit a soft and almost collinear gauge boson from the particle 
I, subject to the infrared cut-off fi on the transverse momentum Q. Note that in 
distinction to a gluon or photon mass regulator, the cut-off /i does not spoil the gauge 
invariance of the theory and can take on arbitrary values, i.e. it is not necessarily taken 
to zero. To logarithmic accuracy, we obtain directly from (0): 

dWfafj?) g 2 s 
alog(/i 2 ) 8ir 2 /r 



The infrared evolution equation (|32|) should be solved with an appropriate initial condi- 
tion. In the case of large scattering angles, if we choose the cut-off to be the large scale 
y/s then clearly there are no Sudakov corrections. The initial condition is therefore 

M{p!, ...,Pn, S) = M B om(Pl, —,Pn), ( 36 ) 

and the solution of ( |3"2"D is thus given by the product of the Born amplitude and the 
Sudakov form factors: 

M(pi,...,Pn;v> 2 ) = M Bo rn(pi, --^Pn) exp (^-^J2Wi(s^ 2 )^j (37) 
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Figure 4: A Feynman diagram yielding Yukawa enhanced logarithmic corrections in 
the on-shell scheme. At higher orders, the subleading corrections are given in factorized 
form according to the non-Abelian generalization of Gribov's theorem as described in 
the text. Corrections from gauge bosons inside the top-loop give only sub-sub leading 
contributions. DL-corrections at two and higher loop order are given by gauge bosons 
coupling to (in principle all) external legs as schematically indicated. 



This is the exponentiation of Sudakov DL in non-Abelian gauge theories [31] . We now 
want to apply this result to the electroweak theory for the subleading Yukawa correc- 
tions at higher orders. Since we are interested here in corrections to order 0(a n L 2n ~ r ), 
each additional loop correction to the universal subleading terms in the previous sec- 
tion must yield two logarithms, i.e. we are considering DL-corrections to the basic 
process like the inner fermion loop in Fig. f|. It is of particular importance that all 
additional gauge bosons must couple to external legs, since otherwise only a subleading 
term of order 0(a n L 2n ~ 2 ) would be generated. All subleading corrections generated 
by the exchange of gauge bosons coupling both to external Goldstone bosons and in- 
ner fermion lines cancel analogously to a mechanism found in Ref. []30| for terms in 



heavy quark production in 77-collisions in a J z = state. Formally this can be un- 
derstood by noting that such terms contain an infrared divergent correction. The sum 
of those terms, however, is given by the Sudakov form factor. Thus any additional 
terms encountered in intermediate steps of the calculation cancel. For the one loop 
process in Fig. |], for instance, we inner include only corrections with top quarks and 
assume on-shell renormalization of the external Goldstone bosons. Thus the correc- 
tions at higher orders factorize with respect to the one loop fermion amplitude and 
A4"Bom"(Pi5 ■■■iPn) — -M ii 00 p (pi > • • • ) Pn) ■ Note that the latter is also independent of 
the cutoff since the fermion mass serves as a natural regulator. In principle we can 
choose the top-quark mass to be much larger than \x for instance. This freedom is 
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not present for subleading terms from gauge bosons, such as the angular contributions 
of the type log | log -^2, which furthermore don't factorize with respect to the Born 



amplitude as mentioned above. Thus, the method suggested in Ref. |32] for the higher 
order angular terms cannot straightforwardly be justified via the non-Abelian gener- 
alization of Gribov's bremsstrahlung theorem 5 . In our case we have for the two loop 
electroweak DL corrections at the weak scale \x = M: 



~.'\v, , ,2, 9 2 ( S ) 



wr(s, m 



1 ' ' Ki~- 



f y 2 
T^Ti + I) + tan 2 6 w ^- 




(3? 



We now want to consider specific processes relevant at future e + e~ colliders and demon- 
strate how to apply the non-Abelian version of Gribov's factorization theorem for the 
higher order corrections. The subleading corrections are then compared to the general 
splitting function approach of section [3[ 

In the case of the amplitude of Fig. f| we must use the quantum numbers of the 
associated Goldstone bosons and we have the following Born amplitude 

e 2 

M Bo m(Pu -iPa) = «7^-(e R l7le^)(/ci - k 2 ) v (39) 

and at one loop we have two fermion loops contributing (ttb and bbt). The renormaliza- 
tion condition is provided by the requirement that the corrections vanish at the weak 
scale, i.e. for s = M 2 , which amounts to subtracting the vertex for that case. The first 
diagram of the two is given by 



■y,Z w 

d n l Tr {uj-fu+df- h)lv{<? + u+ + c^_)(f+ h)} 

(2tt)™ {P -m 2 b + ie)((l + k^ 2 - m 2 t + ie){{l - k 2 ) 2 - m 2 t + ie\ 



+ 5; 



ttb 
ct 



3iQt e 4 m 2 



le^^^ix 8 "-™-^ (40) 

where u± = |(1±75) an d the chiral couplings are given by c£ = Qf for the photon and 

cL = —Qf and c_ = — for Z-bosons respectively. The counterterm 5*i is chosen 

such that the logarithmic corrections vanish for s = M 2 . Thus, the sum of the scalar 
functions is to logarithmic accuracy B23 — = — log-gj. Analogously, we have for 
the bbt quark loop: 



4 2 
p in 

7,^ w 



X 



5 Using the MS-renormalization scheme, however, the subleading pole structure of QCD scattering 
amplitudes at the two loop level is determined only by one loop divergences and renormalization group 
corrections fl3^1. 
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H , k 2 s 

Figure 5: A Feynman diagram yielding Yukawa enhanced logarithmic corrections to 
external longitudinal Z-bosons and Higgs lines in the on-shell scheme. At higher orders, 
the subleading corrections are given in factorized form according to the non-Abelian 
generalization of Gribov's theorem as described in the text. Corrections from gauge 
bosons inside the top-loop give only sub-sub leading contributions. 



d nj Tr {u + (-y)u_(-!f - h)lAc b + u + + c^-)(-f + fe)} 
(2vr)« (P - m 2 t + ie)((Z + h) 2 - m\ + ie)((l - k 2 f - m\ + is) + ct 

:W} "- Ti!) ' •<eR|7le+)(5 23 -^)(fc 1 -fc 2 ) I , (41) 



Adding both results ((5J) and (|4l ) we find 

{(7 2 3 77i 2 5 I 

X ~ 16^2 Jp l0g Jp) (42) 

and the all orders result to subleading accuracy is given by 

M( Pl , ..., Pn ; /i 2 ) = MuooM, -,Pn) exp [~\ £ ^{s, /i 2 )^) (43) 

The subleading Yukawa corrections from the Altarelli-Parisi in Eq. ( p2|) agree with 
the corresponding results from the application of the Gribov-theorem in Eq. (fE|). For 
longitudinal Z-boson and Higgs production, we note that there is only one non-mass 
suppressed elementary vertex with two neutral scalars, namely the Z\H vertex. As 
mentioned above, universal terms are related to the massless limit. For the "Born 
amplitude" of the Higgs- Strahlung vertex we have 

Mlf n =-^-(k»-k») (44) 
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The universal Yukawa corrections to both external \ and H states from an off shell 
Z line are then given by the corrections depicted in the inner fermion loop of Fig. |5|. 
Here we find 



^&i,-,P3) = 3^tx 



UoopU'l) —iPi) AM 2 s 2 

W 

gj Tr { 75 (|Q(f- h)-f + ctuj-KH h)} ZxH 

(2vr)™ (I 2 -m 2 + ie)((Z + fc x ) 2 - mf + ze)((Z - fc 2 ) 2 - m 2 + ie) + ct 

(,V '' ! '^-(^-Bj)^-^) (45) 



16tt 2 s w c w 4M 2 s 2 , 

and thus 

M z u C(Pu -.Pa) = A« {l - ^Y^^log^ \ (46) 



From the same line of reasoning as for the charged Goldstone bosons we find that the 
all orders result is given by Eq. (fj5| ). At the subleading level, this is equivalent to the 
corresponding corrections obtained in Eq. 



5 Top Yukawa corrections for chiral quark produc- 
tion 

In Ref. || all subleading Sudakov logarithms were resummed assuming that all in- 
variants 2pjPi ~ s. The subleading kernel of the infrared evolution equation was 
determined by using the virtual contributions to the splitting functions from QCD and 
applying these results to the high energy regime of the electroweak theory. Soft pho- 
ton corrections were then added by appropriate matching conditions at the weak scale. 
We explicitly restricted ourselves in Ref. to the case where all fermions had masses 
below the weak scale and thus excluded Yukawa enhanced terms. From the arguments 
of section ^ it is now straightforward to include also top- Yukawa terms for chiral quark 
final states. These terms occur for left handed bottom as well as top quark external 
lines. The situation for a typical Drell-Yan process is depicted in Fig. ^| where for the 
inner scattering amplitude we have two contributions. We neglect all terms of order 

1—, — ). Using on-shell renormalization we find for the inner amplitude on the 

left in Fig. ^| for a right handed electron in the initial and a left handed bottom quark 
in the final state from the ± loop for the sum of the 7 and Z contributions: 

4 2 

a a DY ^ 1 4- 1 I — \ 

Aioop = - AsM 2 s 2 d 2 , {eMeR}X 



W W 



_li (h\m-h-k 2 y\h) DY 

(2tt)" (I 2 - ml + ie)((l - h) 2 — M 2 + ze)((Z - k 2 ) 2 — M 2 + ie) ct 
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Figure 6: Feynman diagrams yielding Yukawa enhanced logarithmic corrections to the 
third generation of fermions in the final sate. The inner scattering amplitude is taken 
on the mass shell. No DL-corrections originate from the inner loop. At higher orders, 
the subleading corrections are given in factorized form according to the non-Abelian 
generalization of Gribov's theorem as described in the text. Corrections from gauge 
bosons inside the Goldstone-boson loop give only sub-sub leading contributions. DL- 
corrections at two and higher loop order are given by gauge bosons coupling to (in 
principle all) external legs as schematically indicated. 



4 2 
P 771 

' '£\lA<*)W\h){B*-B%) (47) 



32vr 2 AsM 2 s 2 c 2 

W w 

The scalar functions at high energy evaluate to B 2i — B^ = — log t4j as mentioned 
above. For the diagram on the right in Fig. |6| we have for the bottom again only the 
± contribution. Here we find for the sum of the 7 and Z contributions: 

64DY _ e m tQt / +| I _v 



d n l {hWf\h) 



b c-DY 



1 or 

(2tt)« (P -M 2 + ie)((Z - h) 2 -m 2 + ie)((Z - k 2 ) 2 - m 2 + is) ct 
' (l '" fQt ■{e+\ lv \e^){f L \^\h)(B 23 -B^) (48 



32tt 2 2sM 2 s 2 c 2 



In all cases we renormalize on-shell, i.e. by requiring that the vertex vanishes when 
the momentum transfer equals the masses of the external on-shell lines. All on-shell 
self energy contributions don't contribute in this scheme. For external left handed top 
quarks, the ^ loop is mass suppressed and we only have to consider the x an d H 
corrections. They are given by replacing Q t — > 2Q t {Tf) 2 and Q t — > \Q t in Eq. 



(]48|). It turns out that the Z\H contributions equal the corrections from the 70 and 
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in the case of the bottom calculation. The Born amplitude is given by: 

,2 



M 



DY 
Born 



></l|7"I/r> 



w 



• e 2 

i 9 

set 



l -i^_{ e i\lu\e K ){h\l v \h) ,/l 

\ e ibv\ e K)UK\l U \k) >/R 



(49) 



for top and bottom quarks. In all cases, log ^ terms can be savely neglected to the 
accuracy we are working. Thus we find for left handed quarks of the third generation: 



,Pa) 



- /vt Born 



( g% ^ Tfl 2 S 1 

(Pi. ■••.**) 1 - I^4M^ 6l bg (50) 



For right handed external top quarks we have ± , x an d ^ corrections. In that case 
we observe that the Z\H, 70 ± and Zfi^ loops have an opposite sign relative to the 
left handed case. For the corrections corresponding to the topology shown on the right 
in Fig. ^] we must replace Q t in Eq. (pEBD by Q f — Tj = | for the ± graph. The same 
contribution is obtained by adding the H and \ loops and we find: 



i- 



...,p 4 ) 



- /vt Born 



(Pl,-,P4) { 1 



gr 1 m\ 



16tt 2 2M 2 



M 2 



(51) 



At higher orders we note that the exchange of gauge bosons inside the one loop pro- 
cess is subsubleading and we arrive at the factorized form analogous to the Yukawa 
corrections in section |]. Since these corrections are of universal nature we can drop 
the specific reference to the Drell-Yan process and the application of the generalized 
Gribov-theorem for external fermion lines to all orders yields: 



M(pi 



,;/i 2 ) = A*iioop(pi,...,Pn)exp ^--E w r(s,/i 2 ) 



(52) 



where Wf w 



s, fi 2 ) is given in Eq. fl38|) and the quantum numbers are those of the 
external fermion lines. Since at high energies all fermions can be considered massless 
we can again absorb the chiral top- Yukawa corrections into universal splitting functions 
as in Ref. |J. Thus in the electroweak theory we find to next to leading order the 
corresponding probability for the emission of gauge bosons from chiral fermions subject 
to the cutoff ix: 



s,/i 



9 2 (s) 



16tt 2 



T i (T i + l)+tanX^ log 2 ^ 





'1 + 5 



/.R m f 



M 2 



+ 5 



AM 2 



log- 



(53) 



The second line only contributes for left handed bottom and for top quarks as men- 
tioned above and /' denotes the corresponding isospin partner for left handed fermions. 
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6 Semi-inclusive cross sections 



Up to this point we have only considered the corrections from virtual corrections above 
the weak scale M. The physical photon, however, is massless and must be included in 
a semi-inclusive or fully inclusive way. It is thus necessary to consider now the regime 
for k\ < M 2 . The corrections for external fermion, photon and W ± lines are given in 
Ref. ||, in each case corresponding to the logarithmic probability to emit soft and/or 
collinear particles below the scale M. The high energy solution is then the boundary 
condition for the infrared evolution equation at the scale \i = M. For the longitudinal 
particles, we only have corrections from the charged gauge bosons below the scale M. 
In this regime we also need to consider particle masses. For real photon emission we 
assume that the detector resolution is bounded by /i exp t < M, so that emission from 
real massive gauge bosons does not need to be considered and for simplicity, we restrict 
ourselves here to the soft photon approximation. 

Under these circumstances we are now able to summarize the complete expression 
for observable electroweak cross sections at high energies for all universal leading and 
subleading Sudakov corrections as follows^: 

dcr(j>i, ■ ■ ■ ,p n ,g,g',Hex P ) = da Bo m(pi, ■ ■ ■ ,Pn, g(s), g'(s)) 

( n a n f n 4> ~\ 

x exp - E Wf(s, M 2 ) - £ W/(a, M 2 ) - £ W?(s, M 2 ) 



x exp 



n f 



1=1 



i=l 



- E » 2 ) - m 2 )) - E (<(«, v 2 ) - m 2 



1=1 



1=1 



E^7(m 



i m 3 



t=l 



x exp (w2 xpt {s, rrn, //, /i exp t) 



(54) 



The functions Wf(s, M 2 ) and W/(s, M 2 ) are given in Eqs. ( |23|) and ( |53|) respectively. 
The remaining logarithmic probabilities are given in Ref. and are summarized for 
convenience below: 



W?{s,M 2 ) 




a { S ) n r a '( S ) n/\ , S 

Sw^Po + Sw-^Lft log — 



(55) 



with 



A) 



11 1 

— C 4 n 

12 3 



gert 



24 



9 



-n. 



gen 2A^ h 



(56) 



6 We emphasize that for photon and Z-boson final states the mixing effects have to be included 
correctly as described in Ref. In particular, for transverse degrees of freedom the corrections 

don't factorize with respect to the physical Born amplitude but rather with respect to the amplitudes 
containing the fields in the broken phase. For longitudinally polarized Z-bosons, however, there is no 
mixing with photons and the corrections factorize with respect to the Born amplitude. 
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where n gen denotes the number of fermion generations |34], and rih the number of 
Higgs doublets. Again we note that for external photon and Z-boson states we must 
include the mixing appropriately as discussed in Ref. ||. For the terms entering from 
contributions below the weak scale we have for fermions: 



(4tt) 2 
(4tt) 2 



(log 2 ^ -31og^) , m,«^ 



log 
+ log 



III 

% 

2 a 



l)21og^ 
31og^ 



|i«ffli 



Analogously, for external W-bosons and photons we find: 



(47T) 



log— -l)21og— + log 2 — 



M_ 2 

7z 



M 2 



(57) 



(5f 



1 yflf i M? 

3 ^=1 47r2 iV C 10 g M 2 



E 



3=1 47r 2iv C- 



(59) 



for the virtual corrections and for real photon emission we have in the soft photon 
approximation: 



w expt{ s i m ii A 4 ; Mexpt) 



V n — 2i 
^i=l (4tt 



2 

(4tt)2 

2 

(47r)2 



-log 



2 s 



+ log 2 ^-31o^ 



-"cxpt 



-2 log 



cxpt 



logA-1 21og^ + log 2 



log A - 1 

rn 



, rrii <C 



(60) 



where n is the number of external lines and the upper case applies only to fermions 
since for W we have /i < M. Note that in all contributions from the regime /i < M 
we have kept mass terms inside the logarithms. This approach is valid in the entire 
Standard Model up to terms of order O (log . 



7 Comparison with one loop results 

In this section we compare our results from the infrared evolution equation method 
with the explicit one loop calculation of Ref. J36| for longitudinal scattering in 



e + e~ collisions and the general one loop results from Ref. [M. In Ref. [^] all mass 
singular terms were isolated and the physical basis was used to obtain the DL and SL 
corrections from collinear terms, wave function renormalization and RG contributions. 
The results presented there for fermions (up to Yukawa terms) and transverse degrees 
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of freedom agree with our corresponding results in Ref. 0. Also all terms calculated 
here are, at one loop, in agreement with Ref. JTOf- The results of Ref. 



36 



were 



obtained in terms of the physical fields. We already checked that our method gives 
the correct terms at one loop for transverse degrees of freedom and for fermions (up 
to top- Yukawa terms) in Ref. H. Soft real photon radiation will be included in the 



comparison. This comparison is crucial as mentioned in section £4] since we must check 
that the splitting function approach, in particular the factorization of the DL and SL 

terms takes place with the same electroweak group factor I -^T^T^ + 1) + f^-f- 

from the high effective scalar theory. Only the Yukawa terms factorize differently, 

2 

namely with In the following, the lower index on the cross section indicates the 
helicity of the electron, where el denotes the left handed electron 
the relevant results for ete~ - 



convenience as follows: 

Born r 



da\ 

dn) 




1 + 



W£W£ and etel 



l-2& + 4c 



W^W^ from Ref. 



We summarize 
|H for 



2cls 



2 



2>m\ s s 

2^ l0g i^ + 31 ° g ^ 



+ 2 log 




4 f m l 

Born 




5 - 10cj + 8c 4 , 



65 - Qhcl + 18c 4 



(61) 



l0g ii^ 



s 01 s nl AAE 2 
+ 3 log — + 2 log 



mi 



mi 



+ log 



M 2 



E^log 



i=i 



M 2 



(62) 



The Born cross sections are given by: 





Born 


In, 


L 


da^ 


Born 


dn y 


L 



64tt 2 S 1644 
e 4 1 



sin 2 # 



64vr 2 s 4c 4 



sin 2 



(63) 
(64) 



These expressions demonstrate that the longitudinal cross sections in Eqs. (^) and 
(|5"3] ) are not mass suppressed. Eqs. fl5"T|) and (^) were of course calculated in terms 
of the physical fields of the broken theory and in the on-shell scheme. We denote 



cos6 w and s w = sin6 w respectively. Using e 



v 9 



99 



w 1 



V9 



and 
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= we see that the Born cross section in Eq. ( |63"D is proportional to {g 2 +g 



i2\2 



w 

and Eq. (|63| ) proportional to g' 4 . Below the scale where non-Abelian effects enter, 
we have running coupling corrections only from QED, i.e. g 2 (M 2 ) 



and 



g'\M 2 



where 



'eff 



(M 2 



1 e 2 nf ■ M 2 



(65) 



Thus, the RG-corrections to both cross sections for /x > M are given by (using Ca=2, 



n 



gen 



3 and ri/j = 1 in Eqs. p6|): 
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dtt 



RG 



RG 



da\ 
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da 
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1 + 



z 2 41 



Ucl + 22ci 



^ 1 + 



e 2 41 



644 



log 



M 2 



8tt 2 6c 2 ° g M 2 



(66) 
(67) 



The Sudakov corrections to both cross sections from the infrared evolution equation 
method according to Eq. in the soft photon approximation are given below. The 
quantum numbers are those of the particle-indices and are summarized in Tab. [I]: 
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W ± 


1 





±1 


± 


1/2 


±1 


±1 


X 


1/2 


-1 





H 


1/2 


+1 






Table 1: The quantum numbers of various particles in the electroweak theory. The 
indices indicate the helicity of the electrons and quarks. We neglect all mass terms, i.e. 
consider all particles as chiral eigenstates with well defined total weak isospin (T) and 
weak hypercharge (Y) quantum numbers. In each case, the electric charge Q, measured 
in units of the proton charge, by the Gell-Mann-Nishijima formula Q = T 3 + Y/2. For 
longitudinally polarized gauge bosons, the associated scalar Goldstone bosons describe 
the DL asymptotics. 
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Adding Eqs. (|66"D and (|68|) yields exactly the result in Eq. (|6l|) from Ref. 
Analogously, we have for right handed electrons: 
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(69) 



Again we see that after adding Eqs. (^3) and (|69|) we obtain the result in Eq. (|62l) 
from Ref. [315]. Thus we have demonstrated that to subleading logarithmic accuracy 
our results from the infrared evolution equation method in conjunction with the Gold- 
stone boson equivalence theorem are identical with existing one loop calculations with 
physical fields in the high energy limit. 
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8 Discussion of the results 



In this section we discuss the size of the subleading Sudakov corrections obtained in this 
work. We neglect renormalization group corrections for simplicity and use |- = 

£ = 4jS? = (U3^T28 and € = = ottW The motivation for investigating the 



size of the gauge invariant corrections at the subleading level is two-fold. While this 
discussion is incomplete for processes with a large angular dependence, it is nevertheless 
useful in estimating how good the DL approximation is at higher orders. In addition, 
we gain physical insight into the importance of Yukawa corrections and the partial 
cancellation between subleading terms. 



8.1 Sudakov effects for longitudinal gauge boson and Higgs 
production 

We begin with the Yukawa corrections for external scalars given in Eq. (^3|) with the 
infrared cutoff \i = M. Using the quantum numbers of Tab. [I], we have 
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167T 2 
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1 /I ,\ ? n IV, <? s , , s \ 3 m, 2 , s 



167P 



0.79 log 2 ^ + 4.01 log ^ 



(70) 



where we use M = 80 GeV, m t = 175 GeV and = 0.23. The first thing to 
notice is that the Yukawa enhanced logarithms dominate over the subleading Sudakov 
corrections and enhance the overall Sudakov suppression. At 1 TeV we have log — 
5.05 and thus almost equal contributions from DL and SL terms. At 2 TeV we have 
log — 6.44 and at 3 TeV log -A^ = 7.25. In real calculations, however, one finds that 
the Yukawa terms are always proportional to log-V Since the factor of the Yukawa 
logarithm is uniquely determined by Eq. ( |2"3"D we can replace the respective mass term 
inside the logarithmj] Thus, for (f)^ for instance, we have to consider 
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(71) 



11.24 



With the above mass values we have at one loop about 40 (39) % at 1 (3) TeV from 
the subleading terms relative to the DL corrections and at the two loop level about 
79 (77) % at 1 (3) TeV. The subleading corrections are therefore non-negligible and 

Analogously for \ we can put M = Mz and for H we have M = Mh as arguments of the 
non- Yukawa logarithms in Eq. (123) depending on which mass is the largest in a given process. 
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enhancing the Sudakov suppression. Even at 100 TeV the subleading terms make up 
about 52 % (!) at the two loop level and must be taken into account. The good news 
is that the absolute size of the DL correction per line according to Eq. ( [71]) at the 
one loop level is 5.7 (11.7) % at 1 (3) TeV and at two loops 0.16 (0.7) % at 1 (3) TeV 
relative to the Born cross section. 

The above numbers are valid for both external H and x fields. For ^ we also 
have to consider the purely electromagnetic corrections according to Eq. fl54]) . Thus 
we have on the level of the cross section for each longitudinally polarized IV-boson 
including soft photon radiation: 

-wf(s, /i 2 ) + wY(s, M 2 ) + w2 xpt (s, M, /x, /i expt ) = 

„2 



16tt 2 



1 2 S 

log W 2 



2 log 



I, 2 

f'expt 



(72) 



Thus, the complete size of the corrections for 
chosing // expt = M, is given by 



^ on the level of the cross section, 



■Wf 
9 1 



wj(s, /i 2 ) + wj(s, M 2 ) + <p t ( S , M, fi, M) 



16tt 2 



1.02 log^ + 2.01 log JL 



11.24 



(73) 



It is clear that the DL approximation is much more appropriate for longitudinal W- 
bosons than for the neutral external scalars. For instance we have about 6.5 (4.3) % 
at 1 (3) TeV from the subleading terms relative to the DL contributions. The absolute 
size of the DL corrections relative to the Born cross section at one loop is 7 (15) % at 
1 (3) TeV and at two loops 0.26 (1.1) % at 1 (3) TeV. The subleading terms at the 
two loop level contribute about 13 (8.6) % relative to the DL corrections. 



8.2 Sudakov effects for quarks of the third generation 

In order to estimate the size of the corrections for chiral heavy quark production we 
consider first the case of left handed bottom and top quarks. In this case we have from 
Tab. § 



167T 2 

9 2 



Kh 1 



lVw 2 -i--31og— ^ 



tanX-^log M2 



MV AM 2 



ill 7 



16vr 2 



0.814 log 2 — 7 - 1.246 log — - 1.87 



(74) 



M 2 M 2 

The corrections of purely electromagnetic origin are different for the two cases. In 
general we have 

-w{(s,fi 2 ) +w{(s,M 2 ) +w2 xpt (s } m f} n,M) = 



16tt 2 



, 2 s , M 2 ( , s 
- l0g i^- l0g ^l 2l0g M^ 



2l0g M^ 



(75) 
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The full result for the left handed top quark is therefore given by: 



• (s, M 2 ) - (s, /i 2 ) + (s, M 2 ) + < P t (s, m t , //, M) 



16tt 5 



0.916 log 2 



M 2 



1.769 log 



M 2 



1.39 



(76) 



Thus we see that there is a partial cancellation between the subleading and the Yukawa 
terms and the overall DL suppression is somewhat reduced. In relative terms at one 
loop, the SL corrections are about 33 (20) % at 1 (3) TeV, and at two loops the relative 
size of the subleading terms is 65 (41) % at 1 (3) TeV. The absolute size of the DL 
corrections at one loop is 5.4 (11.1) % per line at 1 (3) TeV. At two loops we have 
corrections of 0.15 (0.62) % at 1 (3) TeV relative to the Born cross section. 

The full result for the left handed bottom quark (with = 4.5 GeV) is given by: 



~Wi(s, M ) — Wi(s, fi ) + Wi(s, M ) + w2 xpt {s, m b} fi, M) 

,2 



16tt 2 



0.829 log 2 — - 1.002 log— -2.31 



(77) 



The partial cancellation between the subleading and the Yukawa terms and the overall 
DL suppression is reduced. In relative terms at one loop, the SL corrections are about 
43 (29) % at 1 (3) TeV, and at two loops the relative size of the subleading terms is 86 
(58) % at 1 (3) TeV. The absolute size of the DL corrections at one loop is 5.9 (12.2) 
% per line at 1 (3) TeV. At two loops we have corrections of 0.18 (0.75) % at 1 (3) 
TeV relative to the Born cross section. 

For a right handed top quark we have from Tab. [I]: 



- W- R (s,M' 2 



167T 2 

9 2 



16vr 2 



0.133 log 2 



3 log 



M 2 



M 2 



+ 



2M 2 
- 3.746 



log- 



m, 



(7f 



Now we need to add the corrections from Eq. (|75|). The full result is thus 

—Wf R (s, M 2 ) - w^(s,fi 2 ) +w^(s,M 2 ) +< pt (s,m t ,/i,M) : 



,2 r 



16tt 2 



!^351og 2 ^ + 1.7371og^ 



3.27 



(79) 



Thus we see that there is a large correction of the top Yukawa terms (a factor of 7.4 
for the relative coefficients) and the overall DL suppression is strongly enhanced. In 
relation to the DL contribution at one loop, the SL corrections are about 69 (60) % at 
1 (3) TeV, and at two loops the relative size of the subleading terms is 139 (120) % at 
1 (3) TeV. The absolute size of the DL corrections at one loop is 1.6 (3.3) % per line 
at 1 (3) TeV. At two loops, however, we have corrections of only 0.013 (0.055) % at 1 
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(3) TeV relative to the Born cross section. Thus, the apparent lack of convergence of 
the DL approximation is irrelevant for practical purposes. 

From a physical point of view it is clear how the large subleading terms can be 
understood. Right handed fermions in general couple only to photons and Z-bosons 
and the coupling to the Z-boson is proportional to (^) for the DL and non- Yukawa SL 

corrections. The Yukawa corrections, however are proportional to {j^j and in addition 
for the right handed top, we have: 

m 2 Y 2 4 

1 - 5.38^ = 5.38- (80) 



2M 2 4 9 

a{M 2 ) ~ 3.35a'(M 2 ) (81) 

The effect is somewhat softened by the electromagnetic corrections from Eq. (|75|). In 
general the size of the subleading terms cannot be neglected at the two loop level for 
all Yukawa enhanced Sudakov corrections discussed in this work. 



9 Conclusions 

In this paper we calculated the universal subleading O (g 2n log 2 " -1 jjp, g' 2n log 2n_1 
logarithmic Sudakov corrections to longitudinal gauge boson and Higgs production to 
all orders. We have employed the infrared evolution equation method and used the 
equivalence theorem to obtain the high energy kernel of the equation for longitudinal 
gauge boson production. All Yukawa enhanced SL Sudakov terms in non-mass sup- 
pressed amplitudes are universal to all orders. This feature is evident in the splitting 
function formalism which we have adopted to calculate the virtual Sudakov correc- 
tions. The approach, concerning in particular the novel Yukawa enhanced subleading 
corrections, has been verified by employing a non-Abelian generalization of Gribov's 
bremsstrahlung theorem. We agree with the literature at the one loop level, which is 
a highly non-trivial check considering the complicated nature of electroweak radiative 
corrections and can serve as an independent confirmation of those results. In addition 
this comparison confirms the validity of the splitting function approach since DL and 
non- Yukawa SL corrections factorize with respect to the same group factor of the effec- 
tive high energy theory. These SL contributions are determined by the spin only and 
are thus identical to those found in a scalar theory with an unbroken SU(2) x £7(1), 
while the Yukawa enhanced SL corrections indicate the spontaneously broken gauge 
symmetry. 

The physical picture which is now emerging is clear: at high energies the SM 
behaves like an unbroken gauge theory up to DL and SL accuracy for fermions and 
transversely polarized gauge bosons. Only Yukawa corrections are novel features in 
this picture. For longitudinally polarized gauge bosons and Higgs scalars, the effective 
theory is given by the Goldstone boson equivalence theorem and contains corrections 
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in analogy to a non-Abelian gauge theory with scalar fields in the fundamental repre- 
sentation. Again, Yukawa terms modify this picture as a unique ingredient of broken 
gauge theories. The mass gap between the electroweak gauge bosons can be included 
in a natural way via the matching conditions in the framework of the infrared evolution 
equation method. Thus all universal Sudakov corrections to DL and SL accuracy are 
known in the electroweak theory to all orders. The remaining corrections which enter 
at this level of precision are given by angular terms of the type log j log -£5. These 
terms are non-universal and don't factorize with respect to the Born amplitude. While 
these terms are known at one loop, for phenomenological applications at future col- 
liders a two loop analysis is desirable. In addition, subleading RG effects of the type 
a n /3olog 2n_1 ^2 coupling effects at higher order should be consistently resummed via 
the inclusion of a running coupling in each loop analogously to the QCD Sudakov form 
factor. 

In summary, all universal Sudakov logarithms in the electroweak SM are known 
at the subleading level to all orders and are non-negligible at future collider energies. 
The inclusion of the leading and full subleading electroweak radiative corrections at 
least at the two loop level will be important in investigating new physics effects at TeV 
energies. 
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